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Abstract. We consider the inverse multiphase Stefan problem, where information on the heat 
flux on the fixed boundary is missing and must be found along with the temperature and free 
boundaries. Optimal control framework is pursued, where boundary heat flux is the control, and 
optimality criteria consists of the minimization of the L 2 -iiorm declination of the trace of the solution 
to the Stefan problem from the temperature measurement on the fixed right boundary. State vector 
solves multiphase Stefan problem in a weak formulation, which is equivalent to Neumann problem 
for the quasilinear parabolic PDE with discontinuous coefficient. Full discretization through finite 
differences is implemented and discrete optimal control problem is introduced. We prove well- 
posedness in Sobolev spaces framework and convergence of discrete optimal control problems to the 
original problem both with respect to cost functional and control. Along the way the convergence 
of the method of finite differences for the weak solution of the multiphase Stefan problem is proved. 
The proof is based on achieving a uniform La^ bound, and LUj^’^-energy estimate for the discrete 
multiphase Stefan problem. 
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1 Description of Main Results 


1.1 Introduction and Motivation 


Consider the general multi-phase Stefan problem ([17]): find the temperature function and phase 
transition boundaries {u{x,t), = 1, J} in £> = {0 < x < .^, 0 < t < T} satisfying the 

following conditions: 

a{u)ut-{k{u)ux)x = f{x,t), {x 

,t) G D, 

u{x,t) y£u^,j = 1, J, 

(1) 



0<t<T, j = l,J, 

(2) 

[w]|x=43(t) = Oj 


0<t<T, j = l,J, 

(3) 

[ k { u ) ux ]\ x=^j (t) = ’ 


0<t<T, j = l,J, 

(4) 

u(x, 0) = <^(x), 


0 < X < £, 

(5) 

k ( u ) Uxlx=0 = 9 { t ), k { u ) ux \ x=i = 


0 < t < T, 

(6) 


where a,k are positive functions on each segment (—oo,m^], [u^ j = 1,..., J — 1 and 
, +oo) with 1st type discontinuity at it = , j = 1, J, where < v? < ■ ■ ■ < are known 

values; each 7 ^, j = 1, J is a known positive number, and [it]|a;=j^ is the saltus of v at defined as 

HU=5, -^1^=?,- 

where (or is the limit value of v at taken in the region {{x,t) : v > u^} 

(or {{x,t) : V < u-’}). In the physical context, / characterizes the density of the sources, (f) is the 
initial temperature, g and p are the heat fluxes on the left and right fixed boundary respectively, 
each represents a phase transition temperature, and (|1]) is the Stefan condition expressing the 
conservation law according to which the free boundary is pushed by the saltus of the heat flux from 
different phases. 

Weak formulation of the multiphase Stefan problem, as well as existence and uniqueness of the weak 
solution to the multiphase Stefan problem was first proved in [25l |33]. We refer to monographs 
miiin] for the extensive list of references. 

Assume now that some of the data is not available, or involves some measurement error. For 
example, suppose that the heat flux, p, at the fixed boundary a; = 0 is not known and must be 
found along with the temperature u and the phase transition boundaries ^j. As compensation for 
not knowing this function, we must have access to additional information, which for instance may 
come as a measurement of the temperature at the fixed boundary x = £: 

u{e,t) = vit), 0<t<T. (7) 

Inverse Multiphase Stefan Problem (IMSP). Find the functions u{x,t), j = 1, J, and 

the boundary heat flux g{t) satisfying ([T])-(l7]). 

The IMSP is not well posed in the sense of Hadamard. That is, if the data is not sufficiently coor¬ 
dinated, there may be no solution. Even if it exists, it might be not unique, and most importantly 
there is no continuous dependence of the solution on the data functions. 

We refer to a recent paper [T] for review of the literature on Inverse Stefan Problems. The one-phase 
inverse Stefan problem (ISP) was first mentioned in [12) . where unknown heat flux is to be deter¬ 
mined under the given free boundary. The variational approach for solving this ill-posed inverse 
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Stefan problem was used in [gno]. The first result on the optimal control of the Stefan problem 
appeared in [?D], where an optimal temperature along the fixed boundary must be determined to 
guarantee that the solutions of the Stefan problem stay close to the measurements taken at the final 
time. In [40], the existence result was proved. In [42], the Frechet derivative was found and the 
convergence of the finite difference scheme was proved, and Tikhonov regularization was suggested. 
Later development of the inverse Stefan problem proceeded in these two directions: Inverse Stefan 
problems with given phase boundaries were considered in [1H HI dU dllll HU [ini 133 [H] ; optimal 
control of Stefan problems, or equivalently inverse problems with unknown phase boundaries were 
investigated in [smiuiMiiisiiiiMiiiiisaiisisiiiBisiiii]. We refer to the monography 
[TH] for a complete list of references of both types of inverse Stefan problems, both for linear and 
quasilinear parabolic equations. 

In two recent papers mm a new variational formulation of the one-phase ISP was developed. 
Optimal control framework was implemented where boundary heat flux and the free boundary are 
components of the control vector and and optimality criteria consists of the minimization of the sum 
of L 2 -norm declinations from the available measurement of the temperature on the fixed boundary 
and available information on the phase transition temperature on the free boundary. This approach 
allows one to tackle situations when the phase transition temperature is not known explicitly, and 
is available through measurement with possible error. It also allows for the development of iterative 
numerical methods of least computational cost due to the fact that for every given control vector, 
the parabolic PDF is solved in a fixed region instead of full free boundary problem. In [I] the 
well-posedness in Sobolev spaces framework and convergence of time-discretized optimal control 
problems is proved. In [5] full discretization was implemented and the convergence of the discrete 
optimal control problems to the original problem both with respect to cost functional and control 
is proved. The main advantage of this method is that numerically at each step, the problem to be 
solved is only a Neumann problem, and not a full free boundary problem. Moreover, the Neumann 
condition replaces the Stefan condition on the free boundary. In a recent paper [3] the Frechet 
differentiability and first order optimality condition in Besov spaces framework is proved and the 
formula for the Frechet gradient is derived. 

This approach is not applicable to multiphase Stefan problem. The reason is that the Stefan condi¬ 
tion on the phase transition boundary includes the flux calculated from both phases. Therefore, it 
can’t be treated as a Neumann condition, even if we include the free boundary as one of the control 
components. In the current paper we develop a new approach based on the weak formulation of 
the multiphase Stefan problem, as a boundary value problem for the nonlinear PDF with discon¬ 
tinuous coefficients in a fixed domain. The main goal of this paper is to solve the IMSP in the 
optimal control framework by employing the weak formulation of the multiphase Stefan problem. 
We prove the existence of the optimal control and convergence of the sequence of discrete optimal 
control problems to the continuous problem both with respect to functional and control. The proof 
is based on the proof of uniform L^o bound, and IF 2 ^’^-energy estimate for the discrete multiphase 
Stefan problem. We address the problem of Frechet differentiability and application of the iterative 
gradient methods in Hilbert spaces in an upcoming paper. 

We describe in Section fL^ the notation of Sobolev spaces which we use in this paper. In Section fTTdl 
we formulate IMSP as an optimal control problem. In Section 11.41 we perform full discretization 
through finite differences and formulate discrete optimal control problem. In Section 11.51 the main 
results are formulated. In Section [2] we prove the existence and uniqueness of the discrete state 
vector. We present the proof of the main results in Section[3] In Section [3d1 we prove Loo estimation 
for the discrete multiphase Stefan problem. In Sectionwe prove IF 2 ^’^-energy estimation for the 
discrete multiphase Stefan problem. Based on these estimations we prove the existence of the 
optimal control in Section 13.31 Proof of the convergence of the discrete optimal control problems 
to continuous optimal control problem is completed in Section [3.41 
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1.2 Notation of Sobolev Spaces 

1 ^ 2 ( 0 , T) - Space of Lebesgue square-integrable functions. It is a Hilbert space with inner product 


r 

{u,v) = / uvdt. 
Jo 


1^00 (0,T) - Space of essentially bounded functions. It is a Banach space with norm 

ll«lkoo[o,T] = esssup |u(t)|. 

0<t<T 

W^iO, T), k = 1,2, ... - Hilbert space of all elements of ^ 2 ( 0 , T) whose weak derivatives up to order 
k exist and belong to L2{0,T). The inner product is defined as 


cT ^ 


{u,v) = 


'0 


s=0 


d^v 
dV dV 


dt. 


L 2 {D) - Hilbert space with inner product 

{u,v) = / uvdxdt. 
JD 


^2 ~ Hilbert space of all elements of L 2 {D) that have a weak derivative in the x direction, 

and such that it belongs to L 2 {D). The inner product is defined as 


/ N f J dudv\ , , 

(u,v) = j^(uv + --)dTdL 


W 2 ’ (-D) - Hilbert space of all elements of L 2 {kl) with weak derivatives of first order, Also 

its weak derivatives must belong to L 2 {D). The inner product is defined as 


, , f / dudv dudv\ , , 


1.3 Multiphase Stefan Optimal Control Problem 

Following the well-known reformulation of the IMSP (see [2I1I33]), we consider the transformation 




v{x,t) = F{u{x,t)) := j k{y)dy. 

ni 

Then = f k(y) dy, = 0 <■■■ < , and our conditions become: 


( 8 ) 


(i{v)vt - Vxx = f{x,t). 

{x,t) e D, 

v{x, t) ^ 

(9) 


0<t<T, 


(10) 

['^] U=Cj(t) = 0) 

0<t<T, 


(11) 

r 11 

0<t<T, 


(12) 





v{x,0) = <^{x) = / k{y)dy, 

0 < X < i, 


(13) 

II 

0 

II 

0<t<T, 


(14) 

II 

II 

0<t<T, 


(15) 

v{t) 

p 




v(£,t) = T{t) = J k{y)dy, 

0<t<T, 


(16) 
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where 


(17) 



and F ^ is an inverse function of F. The function I3(v) is of similar type as a and k. It is positive 


functions on each segment {F{—oo),v^], = 1,..., J — 1 and [u‘^,F(+oo)) with 1st 

type discontinuity at v = , j = 1, J. Now, we can invoke a function b{v) such that b'(v) = P{v). 

Our partial differential equation becomes 



(18) 


Moreover, we’re free to choose the jump of b at the values v = vF We choose them in such a way 


that [6(u)]|^^„i = so that upon integration by parts of (IT^ over D, the integrals over the phase 
transition boundaries cancel out. 

Definition. We say that a measurable function B{x,t,v) is of type SS if 

(a) B{x,t,v) = b(v), V ^ , Vj = 1, J, and 

(b) B{x,t,v) S [b{v^)~,b{v^)~^], v = for some j. 

Note that B{x,t,v) can take different values for different {x,t) when v = for some j. 

Given g, a solution to the Stefan problem (l^- (fT^ is understood in the following sense: 

Definition, v € W 2 '"^{D) fl Loc{D) is called a weak solution of the Stefan problem ffl- dThl) if for 
any two functions B, Bq of type S§, the following integral identity is satisfied: 


T e 


II 


B{x,t,v{x,t))'ift+Vxifx- fif dxdt- 


I 


Bo{x^ 0, ^(x))'i/^(x, 0) dx 


0 0 


0 


T 


T 




0 


0 


Consider the control set 


^R = {g.g& Wl{Q,T), ||5||w,Mo,T] < R}- 


We wish to minimize the cost functional ^ given by 


c/ig) = lk(^,^;ff)-r(t)||i,(o,T) 


( 20 ) 


on where v = v{x,t;g) S Wl’^{D) fl Loc{D) is a weak solution of the Stefan problem in the 
sense of OH). This optimal control problem will be called Problem I. 


1.4 Discrete Optimal Control Problem 


uJr = {4, k = l,n},T = —,tk = kr, uJh = {xi,i = 1, m}, h = —,Xi = ih 


n m 


be grids in the time and space domains, respectively, and we’ll assume from here on that 


m —>■ oo 


as n —>■ oo. 
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Define the Steklov averages 


tk 




Ofc = ^ y a{t)dt, ^ = $(^), (21) 

tfc-l a:; 

tk aii+i 

fik = ^ j J f{x,t)dxdt, k=l,n, i = 0,m-l. 

tk-l Xi 

where a stands for any of the functions p, F, g, or g". Introduce the discretized control set 

n = {[9]nG^-+^:\\[9]um<R} 

where [g]n = igo,gi, ■ ■ ■,gn), and 

n n 

MuWii = 

k=l k=l 

with gkt = ■ Consider now the mappings between the discrete and continuous control sets, 

: W^{0, T) ^ W^{Q, T) as 


1 


tk 


^n{9) = [9]n, for g G ^r, where gk = - / 9(t) dt, fc = l,n, go = g{Q), (22) 


tk-l 


^n{[9\n) = g'^, iOT[g\n€^R] ^"(t) = ^fe-l + (t - tfc_l) , t G [tfc-l, <fe), fc = 1, n. (23) 


Approximate the function h{v) by the infinitely differentiable sequence 


rx+i 


bn(v) = / h{y)uJn{v - y)dy, 


(24) 


where u}„ be a standard mollifier defined as 


/ ^ J Cne |u| < i 

' 0 , L;’* 


(25) 


and the constant C is chosen so that f u}i{u)du = 1. Since b'(v) is piecewise-continuous function 

R 

we also have 


rx+-. 


bniv)= b'{y)ujn{v-y)dy. 


(26) 


Hence is strict monotonically increasing function. Next we define a discrete state vector, which 
represent the solution of the discrete multiphase Stefan problem. 

Discrete State Vector. Given [g]„, the vector function [u([ 5 ]„)]„ = (u(0), u(l), ..., v{n )); v{k) G 
/c = 0,..., n is called a discrete state vector if 


(a) Ui(0) = $i, i = 0,m, 


(b) For arbitrary k = 1,..., n, the vector v{k) G satisfies 

m—1 


'^h {bn{vi{k)))f]i+Vix{k)pix-fikPi - Pfc»7m + = 0, Vr? = (?7i) G (27) 


i =0 
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Given [g]„ the discrete cost functional J^n is defined as 


It 2 

J^n{[g]n) = ^T(vmik) -Tk^ (28) 

fe=l 

where Vm{k) are components of the discrete state vector We define 

■— inf ^ni\g\ri)- 
[sin 65^5 

The discrete optimal control problem will be labeled Problem I„. Furthermore, the following 
interpolations will be considered: 

v{x,t) = Vi{k), x&[xi,Xi+i], i = 0 ,m-l, k = 0,n, 

v{x;k) = Vi{k)+Vixik){x - Xi), x&[xi,Xi+i], f = 0,m-i, 
v'"{x, t) = v{x; k), t G [tk-i,tk], 

v'"{x,t) = v{x;k - 1)+Vt{x;k){t - tk-i), t G [tk-i,tk], k = l,n. (29) 


1.5 Formulation of the Main Results 


Unless otherwise stated, throughout the paper we suppose that 


fGLUD), pGWl{Q,T), <^gW^{0,£) 


(30) 


and in addition that $(a:) is such that the critical values v^,j = 1 ,... ,to are taken by $ on a set 
of measure 0 in the x—space. 

Concerning the behavior of the coefficients a and k at oo, we assume the following assumptions: 


pOO 

/ k{u)du = oo, 

J U\ 

liminf ^ > ao > 0. 
u^oo k[u) 


(31) 

(32) 


Condition dSID guarantees that the domain of h{v) is R. Condition (1321) implies that there is a 
uniform positive lower bound for h' and b'^: 


b'{v), b'„(v) > b > 0, V € R. 


(33) 


for some b > 0. 

Theorem 1. The Problem T has a solution, i.e. the set 


is not empty. 


The proof of Theorem [T] hinges upon showing the weak continuity of the cost functional ^. The 
weak continuity of ^ will be established by proving Loo{D) bound and W 2 ’"^{D) - energy esti¬ 
mation for the solution to the discrete Stefan problem and subsequent use of compactness of the 
family of interpolations. 
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Theorem 2. The sequence of discrete optimal control problems In approximates the optimal control 
problem I with respect to functional, i.e. 


lim Jn, = J*, (34) 

n—>-+cx3 

where 

^ n* — inf ^ ^ 

// [gU G ‘IS chosen such that 


then the sequence g" = t^n{[g]n,,) has a subsequence convergent to some element g^, G weakly in 
14^2^(0,T) and strongly in L2{0,T). Moreover, piecewise linear interpolation ■O’’ of the corresponding 
discrete state vector [t'([g]rae)]n converges to the weak solution v{x,t;g^:) G W 2 '"^{D) of the Stefan 
Problem weakly in Wl'^{D). 

The necessary and sufficient condition for the convergence of discrete optimal control problems to 
the continuous optimal control problem is formulated in [40]. The proof of Theorem [2] is based 
on the proof that the conditions of the general criteria are satisfied. As before, Lqo bound and 
14^2^’^ energy estimation for the solution to the discrete Stefan problem play a significant role in this 
context. 


2 Preliminary Results 


Lemma 1. Given any [g]n G f#", and any h,T, a discrete state vector exists uniquely. 

Proof. First we prove uniqueness. Suppose v and v both are discrete state vectors for a given [g]n. 
Due to (a) from the discrete state vector definition, we have that u(0) = 11(0). For a fixed fc > 1, 
suppose that v{k — 1) = v{k — 1). (l27)) is satisfied for both v and v. Substract the identities for 
rj = v{k) — v{k) to get: 


m—1 

[{bn{Vi{k))i 


2=0 


bn{v^{k))t){v^{k) - 'Oj(fc)) + {vix(k) - Vix{k)y 


= 0 . 


However, 


bn{Vi{k))i - bn{Vi{k))t 


bnjvijk)) - bniVlik - 1 )) 

r 


bnjvijk)) - bn{vi{k - 1 )) 

r 


bnjvijk)) - bn{Vi{k)) 

T 


Thus that the previous summation identity becomes: 


m—1 - 

Y [-{bn{Vi{k)) 

i=0 ^ 


bn{vi{k))){vi(k) - Vi{k)) + {vix(k) - Vix{k)y 


= 0 . 


Since 6„(u) is monotonically increasing, the whole summand is non-negative. Therefore, it is equal 
to 0, which implies that Vi(k) = Vi(k), VI = 0,m. Hence, by induction, v = v. □ 


Now we seek to prove existence. Again we’ll rely on induction. Construct v(0) as given in (a) of the 
Discrete State Vector Definition. Note that ||u(0)|| := max|ui(0)| = max|$i| < ||$||i^roq. Now 

2 i _ ^ 

fix fc > 1, and assume that u(fc — 1) has been constructed successfully so that (l27|) is satisfied for 








all K < k. Moreover, assume that \\v{k — 1)|| < +oo. Notice that the summation identity (l27l) is 
equivalent to solving the following system of non-linear equations: 

= ^b„{vo{k - 1)) + h'^fok - hgji 


-Vi-l 


(vo{k) + ^5„(wo(fc))j - vi{k) = 

(k) + (2vi{k) + - v^+i{k) = ^bnivi{k - 1)) + f,k, i = l,m-l 

-Vm-iik) + Vmik) = hpk 


(35) 


We will construct v{k) by the method of successive approximations. It is critical to remember that 
h, T will be fixed here. Choose = v{k — 1). Having obtained ,we search as a solution of 
the following system: 


-h ^bn{v^+^{k))^ - <(fc) 


!^bn{vo{k - 1)) -h h'^fok - hgl 


-v^_-^{k)+ 2v^^^{k) + ^bn(yf^^{k)) - v^^-^{k) = ^bn{vi{k - 1)) + h'^fik, i = 

-vZ1:l{k) + v^+^{k) = hpk 


(36) 


We now proceed to prove that the sequence converges to the unique solution of (I55|) . Substract 
(I5S1) for N and iV — 1 to get 

Vo^\k) - v^{k) + ^ [bnivo+\k)) - 6„(t>^(fc))] = <(fc) - nf “^(fc) 

2{v^+\k) - uf (/c)) -h ^ [bn{v^+^{k)) - 6„(wf (fc))] = vf^^^ik) - v^_^_-[\k) + vfLi{k) - v^_-[\k) 

v^^Hk) - v^{k) = v^t\{k) - vZ-i(k) 

which is transformed as 

(l + ^C,n) {k)) = <(fc) - v^-^{k) 

(2 + {v^^\k) - vf{k)) = {v(i,{k) - v^_,-,\k)) + {v^_,{k) - vl-,\k)) 


where 




Cn,N--= b’^{9v^+^ik) + {l-e)v^ik))d9, 1 = 0,m-1. 

Jo 


Due to (|33l) . we have Ch n — z = 0, m — 1. Hence we have 




■:r\k) - vr{k) = 

2+“rCn,iV 

- v^{k) = v^t\{k) - v^_^{k) 


(37) 
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Let An := max 5 = (l + I7I') G (0,1). From (1571) it easily follows 

0 < 2 <m \ ^ / 


An ^ SAn—i ^ * * * ^ 

Now it is possible to prove that there exist finite limits 

Vi(k) = lim vf{k), z = 0, 

Af—»+oo 

Indeed, from (1551) it follows that for arbitrary j = we have 

- Ao(5^ < wf+^(fc) - v^ik) < AqS^. 

By summation we have 

+00 M—1 M—1 


(38) 

(39) 

(40) 


+ 00 


- Ao ^ < v^{K) - Ao ^ (5^ < v^{k) < v^{k) + Aq ^ + Aq ^ <5^ (41) 


e=N 


e=N 


e=N 


e=N 


for all M > iV > 0. In particular, by choosing = 0 it follows that the sequence {v^} is bounded 
in Let us now assume 


liminf (k) = lim vf^^’ik), Np < Np+i, p = 0 , 1 ,...; lim Np = + 00 . 

N^+oc p^+oo p—)-+oo 

By choosing in (1411) N = Np we have 


+ 00 




(fc) < uf" (fc) + Ao ^ M> Np 


e=N„ 


which implies that 


+ 00 


limsupz;f^(fc) < r!f^'’(fc) + Ao ^ S^, p=l,2,... 


M—¥oo 




Passing to limit as p —> +00 we have 

lim uf^(fc) < lim wf^’’(fc) = liminf w)^(fc). 

M—¥oc p —>-+oo N^-\-oo 

Since opposite inequality is obvious, it follows that finite limits (1391) exist. 


□ 


Given the existence and uniqueness of the discrete state vector for fixed n, we can uniquely define 
for each fc = 1 ,..., n the vector C,k whose m components Q. are given by 


Ck = [ bn{0vi{k) + {1 — 6)vi{k — l))dd, i = 0,m—l. 

Jo 


(42) 


Lemma 2. Sequence of discrete optimal control problems approximates the continuous 
optimal control problem I if and only if the following conditions are satisfied: 

(1) for arbitrary sufficiently small £ > 0 there exists number Mi = AIi{e) such that cSm(<?) 

for all g S i^R-e and M > Mi; and for any fixed £ > 0 and for all g S '^R-e the following 
inequality is satisfied: 

lim sup (yMi-^Mig)) - ^{g)) < 0. (43) 

M—¥oo ^ ' 
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(2) for arbitrary sufficiently small e > 0 there exists number M 2 = M 2 {s) such that .?^m([5]m) G 
^R+e for all [g]M S and M > M 2 ; and for all [g]M G ^r , M > 1 the following inequality 
is satisfied: 

limsup (^(^m([s']m)) - =^m([ 5 ]m)) <0. (44) 

M—^oo ^ ^ 


(3) the following inequalities are satisfied: 

lim sup (e) > , lim inf (-e) < , 

where ^4±e) = inf ^{g). 

^R±e 


(45) 


Lemma 3. The mappings satisfy the conditions of Lemma\^ 

Proof. Fix e > 0, and let g G ^rs, [g\n = ^n{g)- We observe that 

n 71 f \ 

= 7 / 9{t)dt\ < /5^(t)di = ||5||i,[0.T]> 

^==1 ^==1 V tL ) 0 


(46) 


n ( *fc-i 

^ '^ali = 7 I 7 y 9{t) dt - 5(0) I + ^ 7 7 / 5(i) ^ y 5(i) j < 




/ if: i 


tk-2 

2 


!'«)« 


, CBS 


\0 0 
T t 


^ ^ \tfe-l t-T 

tfc t 

1 


Fubini’s 


< - / / l5(C)l^rfe + E7 / / < 


0 0 




ifc-i i-T 

^'l|2 


- II51Il2[0.t] + Ib'llLsiO.T]- 


From (|lSl) . (IT7l) . we get 

||[5]n|l^. < ll5ll^.(0.T) + II5'IIL(0..) < iR-er + ywlio.ry 
Since g G bF2^(0, T) and r —^ 0 as n —>• 00, we know 

1 ™ llff'llL(o.r) = 0 - 

T3.—^on 1 j 


(47) 


(48) 


Consequently, we can choose r so small that the whole right-hand side of the above equation is 
bounded by R^. By definition then, [g]„ G 
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Now let [g\n G be given and write g" = ^n[[g\n)- We see that 


dg^it) 


dt 


tk 

n p n 

* = H / 9ltdt = '^Tgl^, 

k—1 * k—1 

tk-i 


(49) 


dt = 


y|g”(t)pdt = ^ J (^gk-1 + gkt{t - tk-i)^ 

n n 1 ^ ^ ^ 

= ^ r5^_i + r'^gk-igki + o ^5fc5fc-i + o (50) 




A:=l 


A;=l 




A;=l 


(51) 


Actually, since [g\n S it is the case that ^ where C is a constant independent of n. 

k=l 

This of course implies Tgf,^ < for any k, or equivalently, 

\gk-gk-i\<Cy/T, fc = l,...,n. 

Using m in o , we can write 

T 

P n n ^ n 

/ l5”(i)P dt<Y^ Tgl + Cyff'^T\gk\ + -r^ ^ rgl^ < 






fc=l 


<'^Tgl + CVt^t 

k^l 

Combining (pg|) and ([5^ we have that 

mwi[,,T]<\\[9U\l^,+CVTV^^ 


n - n 

\ fc=i 


T9lt- 


(52) 


fc=i 


\ Z! Z + 0(77). (53) 

\ fc=i 


/c=l 


Owing to (153)) , we can choose n so large that r will be small enough to guarantee that the right-hand 
side will be bounded by (i? + e)^. Hence g^ G ^R+e for all n large enough. □ 


Lemma 4. There is at most one solution to the Stefan problem in the sense of m- 


That a solution to the Stefan problem in the sense of m is unique follows by an argument analogous 
to that presented in Section 9 of Chapter V of m- Indeed, we will prove uniqueness in a wider 
class of solutions than that given in m- Suppose that v G Loo{D) only, not necessarily in the 
Sobolev space h 7’^(I1), and that for any two functions B, Bq of type it satisfies the identity 


T e 


B{x,t,v)ipt-\- vifxx + fil> dxdt + / Bq{x,0,^{x))'iIj{x,0) dx+ 


0 0 


i 1 

Jp(t)'ip{£,t)dt-J g{t)-ipiO,t)dt = 0 , Wip € W 2 ^{D),- 4 ){x,T) = 0 ,V’x( 0 ,t) = tfxi£,t) = 0 . 


(54) 
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The class of functions satisfying the above definition contains the class of solutions given in (IT^ . 
Suppose V and v are two solutions in the sense of (IMl) . Due to our assumption on $, subtracting 
(l54l) in V from that in v guarantees that the second integral in (I54|) vanishes, and we obtain: 


T t. 


{B{x, t, v) - B{x, t, v)) {-ijjt + a{x, t)'ipxx) dxdt = 0 


0 0 


where a(x,t) = — --rr. For (x,t) € D such that v(x,t) = v(x,t), it is the case that 

a{x, t) = 0. Otherwise, since B and B are strictly increasing on v a.e. {x, t) G D, it follows 
that a is non-negative for a.e. {x,t). Moreover, the a.e. positiveness of b'{v{x,t)) implies that 
b = essinf b' > 0 and that b is strictly increasing, and so for almost every {x,t) (assume that 
v{x,t) < v{x,t) for the sake of notational simplicity). 


V — V 

< 

V — V 

v{x,t) 

V 

f b'{w)dw+ (6(u*)+— 5(u*)“) 


J bdv 

v{x,t) i:v'^£{v{x,t),v{x,t)) 




Thus a is essentially bounded, and esssup a{x,t) = ai < -|-oo. Fix e > 0, and take as ip{x,t) the 
solution of the Neumann problem 

iljt +{a{x,t) + = F{x,t), = 0, ip{x,T) = 0, (55) 


where the e is added to ensure the uniform parabolicity of the conjugate diffusion coefficient, and 
F is an arbitrary smooth bounded function in D. Note that (15511 is the conjugate heat equation. 
There exists a unique solution G W^’^^D) of the problem ([55]) ( [57] ). Our goal here is to use the 
arbitrariness of F to obtain that B — B = 0 a.e.; to this end, notice that through the use of (1551) . 
we can write 

T i 


J j [B{x,t,v) — B{x,t,v)) {F — dxdt = Q. 
0 0 


(56) 


Thus our goal will be attained if we have an energy estimate on ijjxx for solutions of (|55l) . In 
the following, we prove a sufficient estimation for the analogous Heat Equation (the result follows 
immediately for the conjugate one by a simple change of variables). Let a^(x,t) = a(x,t) and 
for simplicity we don’t write the superscript. Multiply the non-conjugate version of (1551) by tp^x 
and integrate it over the rectangle Dt := (0,£) x (0,t) to get 


t e 


t i 


t t 


t i 


0 0 


I J~ ^ ~ J j P'4’xxdx dr = j j Fxipx dxdr - J Fijjx 

0 0 0 0 0 0 0 

t t e t e 

J J {{'tpr)x'fpx + aipl^) dx dr - J tprpJx ^dr = J J [{ip ^)rpJx + aipl^) dx dr = J j F^^ip^^dx dr 

0 0 0 0 0 
e i t i t i 

\ j ^l{x,t)dx-^ j ipl{x,0)dx + J J aiPl^dxdT<^ J j {iPl + F^) dx dr 


0 0 


t e 


t e 


0 0 
t i 


J 'tpl{x,t)dx -\-2 J J a'il)^^dxdT< J J dx dr J J F^dxdr. (57) 


0 0 


0 0 


0 0 
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t i 


Letting now = J J ipldxdr, it is clear that y'{t) = J dx, thus that ([Ff)l implies 


0 0 


t £ 


y'{t) < y(t) + J j Fl dxdr. 


0 0 


By Gronwall’s Inequality now (more precisely Lemma 5.5 from Ezl), we deduce from the above 


differential inequality that y{t) < e* — 1 f f dx dr, or in other words 


t i 


0 0 


t e. 


t £ 


ip'^{x, t)dx < e* — 1 


F^ dx dr, Vt G (0, T] 


0 0 


0 0 


t £ 


t I 


(EH 


j ip‘l.{x,t)dx + 2 J J atpl^ dx dr < e* J J F^dxdr, VtG(0,T]. 


0 0 


0 0 


T £ 


The first of the above inequalities implies that ess sup / t)dx < e'^ J f F^ dx dr. Now, since 

0<t<T 0 0 0 

V't = atpxx + F, we have 

II'*/'* llL2(£>t) ~ + ^Iliaf-Dt) — (II ®'*/’a;a: II i2 (-Dt) + ll^lli2(Dt)) — ^ || || (£)j) + 21| T’H 

t £ 

< + 2ao J J dxdr < + aoe*||T’a;|||^(£)^)^ . 

0 0 

These results combined provide the energy estimate we need: 

T £ 

/ / *^2;dt + ess^si^ IIV’:d||L(o.c) < 2(||F||i,(D) + ooe^||F^||i^(£,)). (58) 

0 0 


Having (1551) . we can now observe that 


T £ 


{B — B)e'il;l^ dx dt 


0 0 


T £ 


0 0 


[B- B) -^ (a + e) 2 dx dt 

(a + e)2 


Holder’s 

< 


T £ 


2 / T £ 


< 2esssup b{v) | J J 


^0 0 


(a + e) 


dx dt 


(a + e)(V'L) dxdt] < 


vO 0 


T £ 


< Cy/eesssup b{v) iJJ —^^dxdtj —7> 0 as e —>■ 0 

Vo 0 / 


where C is a constant depending only on T and oq. Recall that e < a + e, and so the integral on the 
right-hand side of the above inequality is bounded above by the area of the rectangle D. Therefore, 
(1561) now implies 

T £ 


[B{x, t, v{x, t)) — B{x, t, v{x, t)))F dx dt = 0. 


0 0 
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Owing to the arbitrariness of F, the above equality implies that B{x,t,v{x,t)) = B{x,t,v{x,t)) 
a..e.{x,t) € D, meaning b{v{x,t)) = b{v(x,t)), a.e. {x,t) s.t. v{x,t) ^ v\j = Since b is 

strictly increasing, we therefore have v{x,t) = v(x,t) a.e. {x,t), so v and v coincide as solutions 
in the sense of (l54l) . and thus we have proven uniqueness in this large class of solutions. □ 

Corollary 1. If a weak solution exists, all of the sets j = 1, J have 2-dimensional measure 0. 

Proof. The proof of uniqueness gives us that Bi{x,t,v{x,t)) — B 2 {x,t,v{x,t)) a.e. on D, for any 
two functions Bi , B 2 of type The functions of type ^ generally differ on the sets , so if one 
of them has positive measure, we arrive at a contradiction to Lemma [H □ 


3 Proofs of the Main Results 


(n) + lbllL„o(o.T) + ll5”lliv|(o,r) + ll‘i’l|Loo(o.f)) (59) 


3.1 Loo-estimation for the Discrete Stefan Problem 

Theorem 3. Suppose that p S S Loo(0,^),/ S Loo{D). For [g]n € and n,m large 

enough, the discrete state vector [f([5]rt)]rt satisfies the following estimate: 

IIHnlU„„ := ^max^ ^mm^ < Coo^H/IU, 

where Coo is a constant independent of n and m. 

Proof. Fix n arbitrarily large. Note maxbi(0)| < 111*1^00(0,^)- Consider a positive function 
7(x) G C^[0,^] satisfying 

7(0) = i, 7(£) = i, y(0) = 1, = i < 7(3;) < 1, X G [0,£]. (60) 


Define 7^ = 7(xi), i = 0,m, and denote as x* the value in [xi,Xi+i] that satisfies (by MVT) 
7(xi+i) — 7(xi) = 7'(x*)/i. Transform the discrete state vector as 


Wi{k) = Vi{k)ji, i = 0,m, k = 0,n. 
System (1551) can be rewritten as: 

' hCk^ot{k) - vox{k) = hfok-Pk 

CkVitik) - Vixx{k) = fik, i = l,...,m-l . 
'^m—l.xik') = Pk 


(61) 


Since 


Vi{k) = —Wi{k), Vit{k) = —w,i{k), 
7* 7i 


Vix{k) =—^Wix{k) + (—] wfik) =—Wix{k) + (— ] Wi+i{k), 
1^+1 \lij X 7 * \l^ 


‘^ixxikf) — WixxifC) -t- 
7»-i 


Wix{k) -b 


iik) = 


1 


''^ixx {k) 


li+l 

1 \ 1 


-) +(- 

LV7*/s V7i/xJ 


w^xik) + — wfik), 

lij XX 


li 




1 \ 1 ^ix T '~)ix 

I — ^ixx “T 7«; 

/ XX 7 j _ i 7 j 7 j+i 
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Thus Wi{0) = 7i$i, 


i = 0,m, and for A: = 1, n, 


^Ck^otik) - j^wo^ik) - {j^)^wo{k) = hfok - 

- ^Wi^^ik) - [(i)^ + (i) J wUk) - {^)^.Mk) = hk, 


i = 1, m — 1. (62) 


Furthermore, transform Wi{k) as: 

Ui{k) = Wi{k)e~^^^ , i = Q,m, k = 0,n 

where 


(63) 


>^ = j{mcio,i] + \\lTc[o,i])^ ( 64 ) 

and if G tfc] satisfies through the MVT that = Xe^*'‘ {tk — tk-i) = Ae'^^^r, then 

Wit{k) = e^*'‘~^Uit(k) + Ae^‘ Ui{k). 


So Wi(0) = Wi{0) = i = 0,m, and for k = l,n, the vector u(k) satisfies the system 

- j^uoAk) + - (i) J uoik) = e-^*^{hfok - gl) 

:^Cfee"^^M*t(fc) - ^^Ui^x{k) - Uix{k) + 


J u,{k) = hke 

^Um-iAk) + ^^m(fc) = e 


(65) 


i = 1, m — 1 




Pk 


Now fix ki < n, and define the following sets of indexes for convenience: 

= {(*, fc)|i = 0,... ,m, k = 0,..ki}, 
o/F = {{i,k)\i = 1,... ,m — 1, k = l,...,fci}, 

A) = {{i^k)\i = 0,k = ^, ■ ■ ■, ki}, 

An = {(F fc)|i = m, /c = 1,..., fci}, 

=!?<) = {(i, A:)|i = 0,..., m, A: = 0}. 

Unless confusion may arise, we omit the subscript to ^ki ■ It is clear that 

U AqU AmD A}). 

If Ui{k) < 0 in AC, then maxMi(A:) < 0. Suppose that 3(i,A:) such that Ui{k) > 0. Then 
max'Ui(A:) > 0. Let {i*,k*) G AC be such that Ui*{k*) = maxMi(A:). 

If (i*,k*) G Ao, then Ui*(k*) = maxy^^i < max$i < max$(x). 

i i [0,d 
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If {i*,k*) € then i* = m, Um-i,x{k*) > 0 and we can choose h small enough that = 

G (—§1—5) so that 


'ymjm-l 


Um{k*) < e 


.m < 


7m7m-l < g-Ai, 


—7'(a;"* 


Pk-- 


If ii*,k*) G %, then i* = 0 ,uot{k*) > 0, uoxik*) < 0. Notice that = -:^7oa;. Note 

lox = so for h small enough, we can ascertain 70a; = 7'(a;°) G (5, f )• It follows 


(h/ofe. 

7o7i 



uo{k*) < e (h/ofe. - 5fc.)- 


If {i*,k*) G then Ui*t{k*) > 0, u^*oox{k*) = -j^{ui»+i{k*) -2uj.(fc*) +Ui*-i{k*)) < 0. For 
{i,k) G the corresponding equation in (1651) is equivalent to 


Ck^ '^ixxik^ 

li li+l 


-) +(- 


.T* \ / XX. 


u^x{k) + 

Ui{k) = fike~^*’°. 


( 66 ) 


Define the sets 




+ - > 0 L ^_ = <^ (i, fc) G ^ 


< 0 


And it’s clear Jk' = jVj^ U ./(A. Suppose {i*,k*) G JVj^. Then owing to (1551) and Ui*x(k*) < 0, we 
can write 


L7z* 



Ui*{k*) < fi-k-e 


— Xtk* 


(67) 


If instead {i* ,k*) G JV-, then we can use dMl) and the fact that Ui*x{k*) > 0 to achieve again 
dSH). Therefore, (|67| is achieved in any case. We can choose r so small that e ^ V/c. 

Observe that 


J_pA(i,.-A*) 

O* 


1 'Ji*x 

~‘^i*xx H ^i*x 


li*+l 


72 ^- 172^+1 


< 


4(2||7"||c[07]) + 16(2||7'|| c[ 07 ]) 


b 1 
64 


< 


Then by (l64ll . it is the case that the coefficient of Ui*{k*) is positive independently of i*,k*. 
Therefore, 




— Xtu 


where is a constant depending only on 7 and b. 
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We can put together the obtained estimations to deduce that for (i, k) G ^ki, 


u,{k) <maxui{k) <ma.x\0, ||«'||c[o.^], ||p||c[o.t], IIs”!! 


C[0,T] 


L^iD), C'7ll/IUcx=(r>) 


}■ 


But because Ui{k) = 'jiC ^^’‘Vi{k), we have the following uniform upper bound for the discrete state 
vector: 


i(fc) < max |o, ll^Hcfo.f], lbllc[o,T]) \\9'^\\cio,T] + \\f\\L^(D), ^■y\\f\\L^(D)'^^ {i,k)€ ^ki¬ 


ln a fully analogous manner, we arrive at a uniform lower bound for the discrete state vector: For 
(t, k^ G ; 

min |o, -||$||c[o_f], -|b||c[o.T]> -|b"llc[o,T] - ll/llL„o(r>)j -C'7ll/lb=o(£')}- 


Combining the uniform upper and lower bounds imply (j59p up to ki. But fci was arbitrary in 
1,..., n. Theorem is proved. □ 


3.2 W 2 ’^- energy estimation for the Discrete Stefan Problem 

Theorem 4. Suppose that p G tF2^(0,T),$ G W2^(0,t),/ G Lao{D). For [g]n G and n,m large 
enough, the discrete state vector [«([(?]„)]„ satisfies the following estimate: 

n m—1 /m—1 \ n m —1 

IIMnlll H + Z! H < (68) 

k—1 i—0 ~ ~ \ 2=0 / k—1 2=0 

< ^oo (ll‘^llw2b0/) ll2'llw2bO,T) + lb llw2b0.T)) 


where Coo is a constant independent of n and m. 


Proof Consider n and m large enough that Theorem |3] is satisfied. In (l27l) . choose ij = 2Tvfik). 
Using (HU), write {bn{vi{k)))i = Qvifik). Also, use the fact that 

2TVixik){vit{k))„; = 2TVix{k){vix{k))i = 

= ViAk) + vLik) - 2viAk)viAk - 1 ) + v^^k - 0 - v^ik - 1 ) = 

= Vixik) + (vix{k) - ViAk - 1)) - v'^Ak - 1) = 

= vfAk) - vA{k - 1 ) + AvAfik). 
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We thus have 


m —1 m —1 m —1 m—1 

S-r Kl(fc) - XI - 1) + X Klt(fc) = 

z—0 i—0 i—0 i—0 

m—1 

= 2r ^ hfikViiik) + 2rpfcx;^i(A:) - 2T^^x;oi(A:). (69) 

i^O 

Estimate the right hand side of (j69l) by applying Cauchy Inequality with e > 0 in the first term. 
Recall that b'^{v) > h, Vu. We will have: 

m — 1 m—1 m — 1 m—1 

S-r X + X Hl(fc) - X - 1) + X Klt(fc) < 

2=^0 z—0 z—0 z—0 

m — 1 ^ m — 1 

< &T X X + ‘^TpkV^iik) - 2Tg^Voi{k). (70) 

i=0 ® i=0 

We can absorb the first term on the right hand side of (l70)l to the first term on the left hand side. 
Hence: 


m — 1 m — 1 m—1 m — 1 

^ X + X - X -1)+X - 

z=0 z=0 z—0 z=0 

^ m — 1 

- X + “^TPkVmtik) - 2Tglvot{k), Vfc = 1, ..., n. (71) 

i=0 

Perform summation of m for k from Itoq, 2 < q < n. The second and third term on the left 
hand side telescope, and we obtain: 

q m — 1 m—1 q m—1 

^X ^ X ^^Itik) + X +X X < 

k—1 z—0 z—0 k—1 i—0 

m — 1 . q m—1 q q 

< X Hl(0) + ^ X X + 2 X ^PkVmtik) - 2 X T9>ot{k). (72) 

z=0 k—1 z=0 k—1 k—1 


Use the summation by parts technique on the p and g sums: 

q q q q q-l 

'^TPkVmtik) = ’^PkVmik) - '^PkVmik - 1 ) = '^PkVm{k) - '^Pk+lVm{k) = 

k—1 k—1 k—1 k—1 k—0 

q-l 

= - '^TpktVmik) +PqVm{q) - PlVmiO) 
fc=l 

'^rg^voiik) = -J^^dktMk) + gqVo{q) - g^voiO). (73) 

fc=l k=l 

In view of (|73|l and borrowing (l59l) from Theorem [3l d72|l yields (through Cauchy Inequality): 
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q m —1 771—1 q m—1 

k—l 1—0 1—0 k—1 i—0 

m — 1 -. q m—1 q—1 q—1 

< X ^ X ^ X + X ^Pkt +X + 

1—0 k—l i—0 k—l k—l 

+ 2i9-l|IHn|l?^ + 2(|b||L^(0,T) + ||ff"'||L<^(0,T))|IHn|koo- (74) 


Through the dehnition of the Steklov average, the Cauchy-Schwarz inequality and Fubini’s Theorem, 
for h small enough we have the following results: 

m— 1 

X ^ II^'I1L(0,^) + W\\hie-h,e) < ll^llw,Ho.^)> 

1^0 

Ilff^llioolO.T) < ll5"llvVi(0,T) - 77) 

= X^(5fc”-)" < ll(5”)'llL(0,r) < llff”frv,M0.T) < 7?^ 

k=l k=2 

9-1 9 

= ^^'''Pkt — \\P IlialO.T) — IbllwjbO.T)) 

fc=l k=2 

q m—1 

'ZrY^hfl^ 

ll/lliaCn)- (75) 

k—l i=0 


Applying the results in (1751) to (1771) . 

q m—1 m — 1 q m — 1 

X ^ X + X +X X < 

k—l 1—0 1—0 k—l 1—0 

< Coo-I--|--I-||g llwjTo.T)) (76) 

where Coo is a constant dependent on b,T,R, but independent of n,m and q. Since, g = l,n is 
arbitrary, from diel), dMl) follows. □ 


Theorem 5. Take set of assumptions ^. Let {[g]n} be a sequence in n such that the sequence 
of interpolations {tf^n{[g]n)} converges weakly to g € 1T2^[0,T]. Then the whole sequence of inter¬ 
polations {D’’} of the associated discrete state vectors converge weakly to v = v{x,t-,g) € 
with V the unique weak solution to the Stefan Problem in the sense of (Gl). 


20 


Proof. By the definitions of the interpolations given in (l2^ . and by using (1^ we deduce that 


II^'’’IIl=o(£') = esssup \v'"{x,t)\ 

{x,t)GD 


max 

0<A;<n 


( 


max |z;z(fe)|) 
0<2<m 7 


lIHralUooJ 


(77) 


T e 




T £ 


0 0 

n — 1 771—1 


dxdt <2 EE Thvf^ik) + ^E E < 2(11 [t 


\n\\g 


2 I ll(j)l|2 


wi(o,e) 


0 0 
T e 


k—0 i—0 

n 771 — 1 


fc=l 1=0 


71 rm—1 


n n I V t t V X - X 

/ / dxdt <2 E E + 3^"^'L-(^)) < 3 E E + 2Thi;^j(fc) 

Q Q fc— 1 1 — 0 fc—1 _ 1—0 

71 71 

'^Thvl^tik) < 2'^Th(vl^_^ t{k) + < 2\\[v]n\\% + 2h^\\p\ 


\WpO,T)- 


(78) 


fe=i 


fc=i 


Since [(;]„ S i#”, then ||5"||wbo,T] < i? + 1 for large enough n. From the energy estimates (l5^ . 
(l68ll and calculations (f^ . (|7^ it is therefore the case that is uniformly bounded in the spaces 
W^^{D) and Loa{D). As such, we may choose a subsequence of that converges weakly in 
1F2^' {D) to some function v G 1^2 ’ [D) fl Loc{D), and thus strongly in L 2 {D), by virtue of which 
we can choose a further subsequence that converges to v pointwise almost everywhere. It is our 
intent to show now that v satisfies the definition of a weak solution to the Stefan Problem. To do 
this, first realize that sequences are equivalent in and sequences {fi} are 

equivalent in L 2 {D), as shown by the following calculations: 


T e 


0 0 


n m-l 

u’’ — dxdt = EE 

k—l i—0 


^tk-l Xi 


71 771 — 1 p T 

= 1 : 1 : / [ 

k—l i—0 


T £ 

EE 

i i k=i i=o 

71 771 — 1 

= E E / - 

k^l 1 = 0 X 


(fi(x; k) — v{x; k — 1) — Vt(x] k)(t — t^-i))^ dxdt = 

T e 

:^(x; fc)(tfc — t)^ dxdt = 7^ y J{vl)^dxdt —^ 0, (79) 

0 0 

{vixik) -Vi^{k- 1) - Vi^t{k){t - tk-i)y dx dt = 


tk-l X 

. tu Xi + i 

71 771 — 1 p p 


^ tk-l Xi 


^ / n 771—1 

t)^ dxdt=-Tij2T^J2 


T £ 



0 0 


71 771 — 1 


Xi + 1 


dx dt = EEM Uk)+vUk)ix-x.) 

1 _1 ^ — n 


k—l i—0 

71 771— 1 


= EE -Th^vf^{k)dxdt —^ 0, 


k—l i—0 


Vi{k) 


2 

dx dt = 


(81) 


as n,m go to + 00 . Accordingly, v'^ ^ v weakly in W 2 ''^{D) and v ^ v strongly in L 2 {D) and 
pointwise a.e. on D along a subsequence. Fix arbitrary if G W2^’^(i7) with if\t=T = 0. Actually, 
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due to density of C^iD) in W 2 ^{D), without loss of generality we can consider ijj G C^{D) and 
= 0. Define V'i(fc) = V* Vfc, and consider the interpolations: 

V(x,t) := ipiik), := := ipit(k), 

Xi<x<Xi+i, tk-i<t<tk, i = 0,m, k = 0,n. (82) 


It is readily checked thatconverge uniformly on _D as n, m —>■ oo to the functions 
respectively. Fix n. For each k in (l27)) as satisfied by the discrete state vector [u([5]„)]„, choose 
rji = Tipiik), Vi and sum all equalities (l27ll over k = 1,..., n. The resulting expression is as follows: 


n m —1 


{bn{v^{k)))-'il;^{k) + v^x{k)'il^ix{k) - f^ki^^{k) \ - + Y^xgkMk)- (83) 


k—1 i—0 k—1 

We transform the hrst term through summation by parts: 

n m—1 n m—1 n m—1 




k—1 i=0 


E-E h{bn{vi{k)))^ipi{k) = EE hbn{vi{k))4>i{k) - EE hbn{v^{k - l))V'j(A:) = 

k—1 i—0 k—1 i—0 

n m—1 n—lm—1 

= EE hbn{vt{k))tf}i{k) - EE hbn{vi{k))ipi{k + 1) = 

k—1 i—0 k—0 i—0 

n—1 m — 1 m — 1 m—1 

= - ^ T ^ hbn{vi{k))ipit{k) + ^ hbn{vi{n))ipi{n) - ^ hbn{vi{0))tpi{l) = 


k—1 2=0 


2 = 0 


2=0 


T-r 


bn{v{x,t))'ljjl dx dt — I bn{^{x))lp'^ {x,t) dx. 

0 0 0 

Thus, (1551) can be rewritten as: 

re e 

— bn{v)^l + v^tl^x — dxdt— / bn{^)'ip^ {x,t) dx— 


0 0 
T 


0 

T e 


— J dt + j {0,t) dt + j J bn{v)il)l dx dt = Q. 


T-t 0 


(84) 


Theorem [3] implies that if := {y G R | 3(a:,t) G D s.t. v{x) = y} (i.e. is the range of v), then 

OO _ 

the set iF = (J is bounded in R, hence its closure is compact in K. Because of the piecewise 

n—1 

continuity of b, it follows that b{v{x, t)) G Loo{D), and therefore ||&n(u)||L„„(D) < C. Since I? is a set 
of finite measure, \\bn{v)\\L 2 (D) < C, so that a subsequence {6„, (w(x, t))} can be constructed so that 
it converges weakly in L 2 {D) to a function b{x,t) G L 2 {D). Through a similar argument, we can 
choose this subsequence so that bm (^(a^)) converges weakly in ^2(0, £] to a function bo{x) G ^2(0, £]. 
Take a diagonal of these subsequences as the whole sequence. We see that 


T i 


T e 


J J bn{v)il}l dx dt < i J J b‘^{v)dxdt\ i J J ii^D^dxdt 


< 


T-r 0 


r—r 0 


— ll^rt('c)l|L2(£>)||'*/'t lli2([0,qx[T-r,T]) 


0 as n —>■ OO. 


(85) 
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Now, due to (1851) . the uniform convergence of ipl respectively to tp, ipx, ipt and weak 

convergence of bn(v), &„($), 5 " to b, Vx, bo, g in the respective L 2 spaces, then as n —)> 00 , 
(l84l) implies 


T i 



0 0 


b{x, t)ipt + Vx'ipx — fpJ dxdt — J bo{x)'ip{x, 0) dx— 

0 

T T 

— J p{t)'ip{£,t) dt + J g{t)ip{0,t) dt = 0. 


It can be checked that both b and bo are functions of type If at the point (x, t) 


v{x, t) —>■ v{x, t) ^ 


( 86 ) 




then we have 

bn{v{x,t)) = j UJi/„{\v{x,t) - u\)b{u) 

v{x,t) — ^ 

On the contrary, if at the point (x, t) we have 

v{x, t) —>■ v{x, t) = 


du —> b{v{x,t)). 


then we have 

b{v^)~ < liminf &„({}(x, t)) < limsupt)) < b{v^)'^. 

n ^■oo n—foo 

Since the sequence {bn{v)} converges to b{x,t) weakly in L 2 {D), by Mazur’s lemma there is a se¬ 
quence of convex combinations of elements of {&n(x)} which converges to &(x, t) strongly in L 2 (D). 
Therefore, there is a subsequence of convex combinations which converges to b(x,t) a.e. in D. It 
easily follows that b = B(x, t, v{x, t)) is a function of type In a very similar way, it is seen that 
bo = i3(x, 0,$(x)) is of type Hence, by definition, x is a weak solution to the Stefan Problem 
in the sense of dm). From Lemma 0] then, v is the unique solution, which implies that v is the only 
weak limit point of the sequence {'C’'}- Therefore, the whole sequence {x’'} converges to v weakly 
in IT 2 ^’^(i:»). □ 


3.3 Existence of the Optimal Control 

Consider a sequence {gi} € such that ^(gi) \ Since {gi} is uniformly bounded in 

IT2^(0,T), it is weakly precompact in Therefore, there exists a subsequence {gi^} which con¬ 
verges weakly in W2^(0,T), say, to g € W2^(0,T) S For ease of notation, take this subsequence 
as the sequence {gi}. 

Let vi = v{x, t; gi) and v = v{x, t; g) be solutions to the Stefan problem in the sense of (fT^ with gi 
and g respectively. Then for fixed I, the sequence of vectors {[5i]n} given by [gi]n = ^n{gi) is such 
that the interpolations g^ = £^n[[gi\n) converge weakly in IF 2 ^( 0 ,T) to gi S W 2 ^( 0 ,T) as n —>■ 00 . 

Therefore, Theorem [5] applies, and so associated to [gi\n the interpolations vj of the discrete state 
vectors converge weakly in 11 ^ 2 ^’ (£>) to vi. As such, 

lhz|lvyi,i(£,) < Ihninf < Clminf (||-f (87) 
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where C is independent of n,m and 1. Thanks to {gi} C ^r, it it is clear from (l5^ and (l68ll that 
the right hand side of (1871) is uniformly bounded. Similarly, one can conclude that 

||w/||L,»(r>) < liminf 

n^oo 

Accordingly, {?;;} £ D L^{D) is a weakly precompact sequence in so that it 

contains a subsequence {vi^} which converges weakly to a function v £ W 2 '^{D)^ and thus strongly 
in L 2 {D). Due to this strong convergence in L 2 {D), a further subsequence of {ujj,} can be extracted 
which converges almost everywhere to D on D. Then the uniform boundedness of this subsequence 
in Lao{D) implies that v £ Loo{D)- Now, take this subsequence as the whole sequence. Each 
of the vi satisfies (HU) with gi and with an arbitrarily fixed function B of type Going to 
infinity along the sequence, we have that we can replace gi with g and vi with v in (HU). Indeed, 
B{x,t,vi{x,t)) — >■ B{x,t,v{x,t)) a.e. on D because of Corollary [T] and the fact that vi ^ v a.e. on 
D. Consequently, D is a solution to the Stefan problem with g. But, due to uniqueness of such a 
solution, it follows that ?; = {) in W 2 '^{D) D Loo{D)- Next, note: 


lim 

l—^OO 


/{g) - c/{gi) 


lim 

l—¥00 

lim 

1—^00 

lim 

1—^00 


lim 

l—^OO 


\\v{£,t) r( 0 llz, 2 [ 0 .T] r(f)||^2[0.T] 

< u - r, -u - r >L2lO,T] - <Vi-T,Vi-T >L2lOX 

\\v{e,t) - Wz(^,i)llL[0.T] + 2 < -u - - r >L2lO,T] 


J \v{£,t) - dt + 2 J - vi{£,t))[vi{£,t)-T{t)) dt 

0 0 


( 88 ) 


By the weak convergence of the sequence {?;;} to v in W 2 '^{D), it follows that we have strong 
convergence in the space of traces. In particular, the integrals in (j88|) vanish as Z —>■ oo. Hence 
lim ^ {gi) = ^{g)- This limit is unique though, therefore it is the case that ^ {g) = so that 


g £ 1^*. 


□ 


3.4 Proof of the Convergence of Discrete Optimal Control Problem 

The proof of Theorem [5] is split into three separate lemmas, as shown below. 

Lemma A. Let = inf ^{g), e > 0. Then lim = lim #*(—£). 

^ R±e £—>-0 ^ £—>-0 ^ 

Proof. The proof of this lemma is very similar to the analogous lemma from [1]. If 0 < < ^ 2 , 

then 

^*(£2) < ^*(ei) < J*< ^*(-£1) < ^*(-£2) 

Hence lim and lim e) > Choose g^ £ ^R+e such that 


lim [^/{ge) - = 0. 

Since {ge} is weakly pre-compact in W2^[0,T], there exists a subsequence e' such that g^' — >■ g* 
weakly in Hd)[0,T1 as e' —>■ 0. Since / is weakly continuous, ^iqe') —t ^iq*) as £ —>■ 0. Hence 
^ ^{g,) > as £' ^ 0. Thus lim 
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From the other side, by Theorem [T] we know G s.t. ^{g*) = If g* G ‘^R\d^R, then 
de* > 0 such that g* G ^rs, Ve < e*, and in this case e) = Ve < e*. If g* G d^R^ 

then 3{ge} with g^ G '^R-e such that g^ —>• g* in W 2 ^[ 0 ,T] as e —>■ 0. The continuity of ^ gives 
us that lim ^{ge) = ^{g*) = Since on the other hand, ^{ge) > ^*(—e), it follows that 


Lemma B. 'ig G '^r, lim Jn{-^n{g)) = /{g)- 


Proof. Take g € ^r arbitrarily. If ^{g) = [g\m and g” = ^n([ff]n), then g"^ ^ g strongly in 
IF 2 ^( 0 ,T) as n —>■ oo. Applying Theorem [5l we have that the interpolations v'^ of the discrete state 
vectors [i^([ 5 ]n)]n converge to z; = v{x,t;g) weakly in IF 2 ^’^(D) as n —>■ oo, and thus the traces 
converge strongly in L2{0,T) to trace By calculations (17^ and (|80)) . the sequences 

{v'^}, {D'^} are equivalent in W^'^iD), so that the v'^{£,-) traces too converge to v{£,-) strongly in 
1 ^ 2 ( 0 , T). If we define 


I 


ik 


T(t)=Tk = - / T{t)dt, tk-i<t<tk, fc = l,n, 


(89) 


— 1 


then f —>■ r in L2{0,T) as n —>■ 00 . Therefore, 


l^ui^nig))-/{g)\ = 


n p 

^T(um(fc) - Tfc)^ - j (v{£,f)-T{f)Y dt 

k=i { 

< ■) ~ OllialO.T] + l|r ~ y] + 


+ 2 


( 1 ^ establishes the lemma. 


yi£,t) - v(£,t)\\v{£,t) - f(t)| + |r(t) - f{t)\\v(£,t) - r(<)| 


0 as n —>■ 00 . 


dt 

(90) 

□ 


Lemma C. For an arbitrary sequence {[g]n} such that [g]n G 

{^n{[g]n)) - yn(,[g]n)'j = 0 . 


Proof. Let 5 " = l3^n{[g]n)- This sequence is weakly precompact, so that a subsequence g"' 
converges weakly to a function g in W 2 ^( 0 ,T). Take this subsequence as the whole sequence. Note 
that 


1 ^( 5 ”) - ^n{[g]n)\ < 1 ^( 5 ") - ^{g)\ + \^{g) - ^n{[g]u)\. (9i) 


Since ^ is weakly continuous, {g'^) — ^{g)\ —>■ 0 as n —>■ 00 . It remains to show that the second 
term on the right hand side of m goes to 0 as n —>■ 00 . Actually, the proof of this fact flows in a 
manner very similar to the proof of Lemma B. From (lOTl) it follows that lim ( ^i£S^ni{[g]ni)) — 

ni—¥oo \‘^ 
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’^ni{[g]ni)j = 0. However, the subsequence chosen was arbitrary. Therefore, the same result 

can be achieved for any subsequence {gn^} of {ffn}- It is then the case that the whole sequence 
i^ni[g]n)) - <^ni[g]n) converges to 0 as n -)> 00 . □ 
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